Scaling of energy spreading in strongly nonlinear disordered lattices 
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To characterize a destruction of Anderson localization by nonlinearity, we study the spreading 
behavior of initially localized states in disordered, strongly nonlinear lattices. Due to chaotic non- 
linear interaction of localized linear or nonlinear modes, energy spreads nearly subdiffusively. Based 
on a phenomenological description by virtue of a nonlinear diffusion equation we establish a one- 
parameter scaling relation between the velocity of spreading and the density, which is confirmed 
numerically. From this scaling it follows that for very low densities the spreading slows down com- 
pared to the pure power law. 
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In linear, disordered, one-dimensional lattices, all 
eigenmodes are exponentially localized due to Anderson 
localization [1[. These models have been first suggested 
for disordered electronic systems |2j, but they are also 
applicable to a wide range of wave phenomena (optical, 
acoustical, etc.) in disordered media [3(. Localization 
effectively stops spreading of energy in such situations. 

By considering waves of large amplitudes, one faces 
nonlinearity and meets the question whether the local- 
ization is destroyed due to a nonlinear interaction of 
eigenmodes. Although this question has been addressed 
numerically [4H2J], mathematically [8|, and even experi- 
mentally in BECs [9[ and optical waveguides [lOj, a full 
understanding is still elusive. It is easier to understand 
how nonlinearity destroys localization leading to ther- 
malization and self-transparency in short random lat- 
tices [11], than to analyze asymptotic regimes at large 
times in long lattices. The most striking effect observed 
in numerical studies is a subdiffusive power-law spreading 
of energy, lasting until maximally available integration 
times. Whether this spreading is an asymptotic state, 
or transforms into a much weaker logarithmic spreading 
suggested by some theoretical estimates [g, [l9( , or even 
stops, remains a challenging problem. 

In this letter we attack the underlying spreading mech- 
anisms by utilizing scaling arguments. The concept 
of scaling has been extremely successful in the under- 
standing of Anderson localization 



121, as well 



descriptions of noncquilibrium phenomena like surface 
growth [l3|. In this paper we demonstrate that the 



spreading of energy in strongly nonlinear disordered lat- 
tices satisfies scaling relations. 

The starting point is a heuristic description of the 
spreading by the nonlinear diffusion equation (NDE) |7( : 



dp 

~dt 



dx V dx 



with 



pdx = E. (1) 



Here E is the total energy (which is conserved) and p 
represents the energy density. Heuristically, the NDE 
describes diffusion that appears solely due to nonlinear- 
ity, for the nonlinear disordered lattices this can be at- 



tributed to a random exchange of energy between modes 
due to their chaotisation. The NDE has a self-similar so- 



lution [14[ describing asymptotic subdiffusive spreading 
with the edge of excitation propagating according to: 



X 
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(2) 



where to accounts a time shift depending on the pec- 
ularities of the initial state. At the moment, it is not 
possible to derive the NDE for a particular nonlinear dis- 
ordered lattice, but one can check if the scaling predicted 
by Eq. (|2|) holds. In other words, one can check if a par- 
ticular nonlinear lattice belongs to a NDE-universality 
class. The main property of this class is that dependen- 
cies on the total energy E and on time are described by 
one parameter a. This allows us to write the law of edge 
propagation in a scaled form 



X 

~E 



t-tg 
E 2 



l/(2+a) 



(3) 



Still, this expression contains an unknown, non-universal 
constant to- We get rid of it by considering the local 
inverse velocity of the edge dt/dX for which holds 



1 dt 
XdA 7 



(4) 



We notice that on the r.h.s. the global density of the 
field w = E I X appears. This allows us to rewrite and 
generalize ((4]) in the form of an one-parameter scaling 
relation 



a(w) = — 



dlogu; 



(5) 



where a dependence of the index a on the global den- 
sity w would indicate deviations from the pure power- 
law scaling given in eq. ([3]). Relations ([2I4I5P define the 
scaling laws to be checked for particular systems; if they 
are satisfied, then we say that the system belongs to a 
NDE-universality class. 



We apply the scaling relations (|2I4I5P to Hamiltonian 
lattices with strongly nonlinear coupling: 
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(6) 



with k > 2, A > 2 being the powers of the local and 
the coupling potential, respectively. Random frequencies 
ui k of the fc-th oscillator are chosen to be independent 
random numbers uniformly distributed on the interval 
(0, 1). Such lattices without local potential (W = 0) pos- 
sess traveling compact waves (compactons) [l5||. With 
a regular local potential the compactons are long-living 
objects dominating the energy spreading and making 
it nearly ballistic, while the presence of disorder sur- 
presses them and makes spreading subdiffusive as will 
be presented below. Parameter W describes the disor- 
der strength whereas /3 governs the nonlinear coupling. 
However, in the case k ^ A by applying a transfor- 
mation q k ->• W a /3- a q k , p k ->• W Xa / 2 p- Ka / 2 p k and 
H -> W Xa l3~ Ka H, with a = 1/(A - «), we can set 
parameters W and /3 to one. Thus, the only relevant pa- 
rameter is the total energy E. Hence, for k ^ A, varying 
the disorder strength W or the nonlinear strength j3 is 
strictly equivalent to appropriate changes of the energy E 
in this system. 

For the special case k = X the system exhibits full 
energy scaling which is seen from the invariance of the 
equations under the following transformations: 



(7) 



Furthermore, we can replace q k —> W~ 1 ' K q k and 
j3 — > j3/W and finally end up at the Hamiltonian © with 
W = 1 and X = k. Now the only parameter is j3, which 
describes the relative strength of the local and the cou- 
pling potentials. 

The case k = X is special because here we can establish 
an exact relation between the parameter of NDE a and 
the nonlincarity index k. Indeed, excluding 7 from the 

2—K. 

scalings of time and energy in (J7J) we obtain t ~ £~ . 
Comparing this with the scaling t — to ~ X 2+a E~ a that 
follows from Eq. ([2|), we find a = 3=p (Note, that we 
have not derived the NDE from the Hamiltonian, but 
solely use scaling arguments to find the exact correspon- 
dance between a and K. Validity of the approach has still 
to be checked numerically below.) From the expression 
for a, we find the spreading law: 



q^q' 


= 79) 


p^p' 


= 7 «/a„ 


t -> t' 


= t/r /2 -\ 


E^E' 


= 7 K £. 



X~(t-to)*?- 



(8) 



For this lattice with homogeneous nonlinearity (in the 
sense that local and coupling potenials have the same 
nonlinearity index k) the one parameter scaling predic- 
tions (|3I4I5|I arc trivially fulfilled with index a being inde- 
pendent on density w, as the energy dependence follows 
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FIG. 1. Energy profiles E k for k — 2 and A = 4 of initially 
localized states at different times t = 10 4 , 10 6 , 10 8 (inner to 
outer, one specific disorder realization, E — 1). 

from exact rescaling. For lattices with nonhomogeneous 
nonlincarity A ^ k the density dependence (|3I4I5|) is non- 
trivial. 

For the Hamiltonian ©, defining the local energy at 
site k as 

Ek = V l + Wulf + ^ ((flfc+x - q k ) X + (q k 9fc -i) A ) , 

we can interprete it as the distribution p(k, t) = Ek, time 
evolution of which we compare with predictions of NDE. 
The numerically obtained profiles of E k (t) are depicted 
in Fig. [T] The peculiar property of strongly nonlinear 
lattices of type © is that the field has very sharp edges: 
one can estimate that the tail decays supcrcxponentially 
(like for compactons in systems without disorder [151]). 
This corresponds well to the property of the self-similar 
solution of the NDE to have a sharp edge, and allows 
us to compare predictions of the scaling theory with the 
numerics for the lattices. 

We start with testing the approach by applying exten- 
sive numerical simulations to the case k — X — 4 where 
theory (JSJ) predicts a = 1/4. Wc integrated the equations 
of motion by means of a fourth order symplectic Rungc- 
Kutta method [16|. In all simulations in this paper we 
used time steps in the range 0.01 — 0.1 that assured 
conservation of energy with accuracy ~ 10 -9 , presented 
results are averages over hundreds of realizations of dis- 
order. Our main quantity of interest is the time AT(L) 
required to excite the next oscillator with L already bee- 
ing excited. This quantity can then be interpreted as the 
inverse velocity of the edge that enters Eq. (|4|5[) . where 
the transition L — > L + 1 implies dX = 1. We defined an 
oscillator to be excited if E k > 10 -50 . (Due to the sharp 
edges of the states this is a reasonable, though arbitrary, 
value; changing this threshold to, say, 10 -20 or 10 -100 
produces similar results.) To obtain the mean value of 
AT(L), we averaged, for each /3, log(AT) for each L over 
disorder realizations. For a better visual display of the 
excitation times on a logarithmic scale we also averaged 
AT over neighbouring sites L. The results are shown in 
Fig. [2] (a), at large L they are in perfect correspondancc 
with the theoretical prediction AT ~ L 5 / 4 . 

We compare these results with a more traditional ap- 
proach, where the width of the field distribution is av- 
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FIG. 2. (color online) Excitation times AT(L) (a) and par- 
ticipation number P(t — to) (b) of spreading states in the 
hamiltonian lattice ([6} with k = A = 4. Values to are ad- 
justed to have a maximally extendended range of power law 
behavior. The dotted lines correspond to the scaling results 
O = (k - 2) /(2k) = 1/4 and have slopes 5/4 for AT and 4/9 
for P. The inset in (a) shows the instantanious time depen- 
dent exponent a from AT ~ L a saturating at the expected 
value a = 5/4. The inset in (b) shows the rescaled spreading 
due to the scaling expectation P/(t — to) vs. t — to- 



cragcd over realizations of disorder for fixed moments of 
time. We have calculated two measures of the width - the 
squared mean displacement (Afc) 2 and the participation 
number of the energy distribution P, as already used in 
literature [a, LJ], [l7|, lllf . However, both quantities behave 
identically in our numerical simulations, hence we show 
only the participation number exemplarily representing 
the spatial extent, defined as P _1 = ^2 k (Ek/E) 2 . Obvi- 
ously, we have P ~ Afc ~ L ~ X. Values of P have been 
averaged over increasing time windows and disorder real- 
izations, the results are shown in Fig. [3 (b). These time 
evolutions fit nicely Eq. ([3j with the theoretically pre- 
dicted value a = 1/4. Hence, we have found the NDE to 
reproduce the correct spreading behavior of this nonlin- 
ear disordered lattice. Based on the assumption of the 
validity of the NDE we were able to calculate the correct 
spreading exponent analytically. Note, that no further 
assumptions or parameter fits were required in this case. 
Although the two methods used - propagation times 
AT(L) and mean packet widths P(t) (or, cquivalcntly, 
L(i)) - appear to be nearly equivalent, the former one 
has two clear advantages. First, it does not possess an 
arbitrary parameter to as the time differences are calcu- 
lated. The second, and more important, advantage is 
that by fixing the system length L in the averaging over 
disorder we in fact fix the characteristic energy density w. 
On the contrary, by averaging the width at a certain time 
we do not fix the energy density, as at a given time the 
variations of the density in different realizations of disor- 
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FIG. 3. (color online) Panel (a): Excitation times AT of 
spreading states for several energies E — 0.05 ... 2. The re- 
sults shown in this plot are averaged over logarithm inter- 
vals on L (see text). Panel (b) shows the rescaled quanti- 
ties log 10 ( AT/L) vs. — log 10 w with w — E/L according to 
Eq. Q; here the solid line is the slope a(w) of the fitting 
parabolic curve according to Eq. l[5|. 

der may be enormous. For the model with homeogeneous 
nonlinearity n = X this is not essential as the time scales 
with energy in a trivial manner. For the nonhomogeneous 
Hamiltonians, to be considered below, this is crucial, as 
one expects the properties of the spreading to depend in- 
trinsically on the energy density, but not explicitely on 
time. 

Above, we have checked the approach on the homoge- 
neous Hamiltonian model with k = A, where the scaling 
with the energy is trivial. Now we apply our method to 
the mostly nontrivial nonhomogeneous case k^ A, where 
the theory based on the NDE predicts one-parameter 
scaling laws (|3|4|5jl . More precisely, we focus on the case 
k = 2 and A = 4 which resembles the widely studied prob- 
lem of the discrete Anderson nonlinear Schrodingcr equa- 
tion (DANSE): its Hamiltonian in the cigenmode repre- 
sentation also possesses a quadratic local disorder term 
and a nonlinear fourth order mode-to-mode coupling. 

First, we investigated the excitation times AT(L) - 
the results are shown in Fig. [3] (a). After performing the 
scaling according to ([4]), all the curves collapse to one 
as seen in panel (b) of Fig. [3] The same approach ap- 
plied to the participation number also leads to a collapse 
of data when the scaling representation (|3|) is used (see 
Fig.0]). The collapse of data for different energies proves 
numerically that the one parameter scaling suggested by 
Eqs. (|4|5|) works nicely for the strongly nonlinear lat- 
tice ([6]) . The validity of this scaling means that asymp- 
totically the spreading in such systems is governed solely 
by the average energy density w. This was assumed in 
most of the previous works on this topic [5|, [D, UM > but we 
present here the first direct numerical evidence of this. 

Fig. [3] (b) shows that within the studied range of two 
decades of variations of the density, parameter a(w) (the 
slope of the curve in Fig. [3]) is not a constant, but a 
growing function of inverse density. In particular, data 
in Fig. [3j can be well fitted with a linear dependence 



a(w) 



-0.3 — 1.51og 10 iy 



(9) 
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FIG. 4. (color online) (a): Participation Number P of initially 
localized states for several different energies E — 0.05 ... 4 
(cf. fig. [3]). Panel (b) shows the rescaled quantities P/E 
vs. (t — to)/E 2 according to eq. (J3J). The red (grey) line 
in (b) visualizes the analytic expectation P/E ~ X/E = 1/w 
obtained from Eq. (0} and (|9]). 



While wc studied in details the strongly nonlinear 
Hamiltonian lattices, it remains a challenge to extend 
the results to lattices with linear coupling terms, e.g. on 
the nonlinear disordered Schrodingcr lattice. The main 
issue here is that for latter situations one cannot define a 
sharp edge of the spreading wave packet, thus the calcula- 
tion of the edge velocity, entering the scaling relation ([5]) 
is problematic. Further studies on lattices of nonlinear 
oscillators coupled by a nonlinear ity of different order, 
e.g. k = 4, A = 6, are currently pursued. 

We thank S. Flach, D. Shepelyansky, S. Fishman and 
S. Aubry for useful discussions. Support from DFG 
project PI220/12 is acknowledged. 



what corresponds to a parabolic fit for the dependence of 
log 10 AT/ L on log 10 w. This means that the spreading of 
energy in the lattice is not a pure power law, but slows 
down as the density decreases. In the intermediate range 
of densities the parameter is close to a = 3, what means 
that in this range the width of the wave packet spreads 
as X ~ t 1 ' 5 , i.e. with approximately the same index as 
found numerically for the nonlinear Schrodingcr lattice 
with disorder in Refs. |5H7(. However, for this model 
no numerical slowing down of the spreading have been 
reported, although a recent theoretical estimation in [19[ 
gives a sub-power law asymptotics logP ~ log ' t. Wc 
stress that an application of the one parameter scaling 
(JSJ) with the empirical law ([9]) also demonstrates a good 
agreement for the participation number as seen in Fig. [4] 
(compare the line in panel (b) with the markers). 

In conclusion, we have studied subdiffusive spreading 
of energy in Hamiltonian lattices with both linear and 
nonlinear disorder and nonlinear nearest neighbour cou- 
pling. Our main result is the one-parameter scaling rela- 
tion, Eq. ([!]), which relates the average velocity of spread- 
ing with the energy density. This scaling relation is mo- 
tivated by using the nonlinear diffusion equation as a 
phenomenological model for the macroscopic properties 
of spreading. We studied in details two sets of exponents 
on the Hamiltonian lattice, k = A = 4 and k = 2 , A = 4. 
In the first case there is no nontrivial dependence on the 
energy, what allowed us to find the spreading index an- 
alytically and to confirm it numerically. This agreement 
of analytical and numerical results is not surprising, but 
still remarkable because it shows that indeed the NDE 
is an appropriate framework to approach systems with 
disorder and nonlinearity. The latter case of linear disor- 
der and nonlinear coupling is mostly nontrivial, here our 
approach gave a density dependent index a(w), that in a 
large range of densities is close toa~3 but grows as den- 
sity decreases in course of spreading. This dependence 
a(w) has not been observed before (e.g. in the DANSE 
model), and it is a first indication of a deviation of the 
spreading from the perfect subdiffusive power law. 
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